Exact ground states of a frustrated 2D magnet: deconfined fractional excitations at a first order 

quantum phase transition 
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We introduce a frustrated spin 1 /2 Hamiltonian whicii is an extension of tlie two dimensional J\ — J2 Heisen- 
berg model. The ground states of this model are exactly obtained at a first order quantum phase transition 
between two regions with different valence bond solid order parameters. At this point, the low energy excita- 
tions are deconfined spinous and spin-charge separation occurs under doping in the limit of low concentration of 
holes. In addition, this point is characterized by the proliferation of topological defects that signal the emergence 
of Z2 gauge symmetry. 

PACS numbers: 71.27.+a, 71.28.+d, 77.80.-e 
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Frustrated magnets are the focus of considerable attention 
because exotic quantum effects are expected to emerge from 
the competition between two or more opposite tendencies. 
While several models in this category are solvable in one di- 
mension, the list is much smaller for higher dimensions. One 
of the most studied frustrated magnets is the spin 1/2 Heisen- 
berg model with first and second nearest neighbor interactions 
Ji and J2. In one dimension, this model exhibits a quan- 
tum transition as a function of Jij J\ from a critical state with 
quasi-long range antiferromagnetic (AF) order to a dimerized 
phase. Moreover, the exact dimerized ground state has been 
obtained for the point J2I J\ = 0.5 by Majumdar and Ghosh 
|[ll]. In contrast, two dimensional (2D) frustrated magnets like 
the Ji — J2 Heisenberg model on a square lattice still hold 
many secrets. Different approaches predict a transition be- 
tween a Neel ordered state and a gapped (non-magnetic) quan- 
tum phase for the region 0.4 < J2/J1 < 0.6. However, the 
nature of this phase is still debated. More precisely, the ques- 
tion is whether it is a uniform spin liquid IH or a spatially 
ordered valence bond crystal j^jsllaSl- 

The interest in frustration and magnetism is motivated not 
only by this widely debated question. There are reasons to be- 
lieve that frustrated magnets may exhibit fractionalized exci- 
tations similar to those which appear in the fractional quantum 
Hall effect. The interest in phases with fractionalized excita- 
tions is generated by the increasing number of experimental 
results showing new physical behaviors in strongly correlated 
systems. For instance, the normal state of the high tempera- 
ture superconductors does not exhibit electron-like quasipar- 
ticles in its spectrum. Recently, different scenarios were pro- 
posed for the realization of points with deconfined fractional 
excitations. Senthil et at |8| proposed that deconfined crit- 
ical points may describe quantum phase transitions between 
Neel ordered states and valence bond solids (VBS) in frus- 
trated magnets. In a second scenario based on the study of 
models for quantum dimers, the deconfinement point sepa- 
rates two VBS and the spectrum of fractional excitations con- 
sists of spinless particles ("photons") with a quadratic disper- 
sion Very recently, Tsvelik JlOII argued that none of the 
previous scenarios is realized for a family of frustrated spin 
Hamiltonians called confederate flag models As we will 
see below, his alternative scenario for the deconfinement point 
at the transition between two VBS has many analogies with 



the case that we study in the present paper 

In addition, during the last decade attention has focused on 
the study of inhomogeneous structures that are proposed to 
emerge from competing interactions. These textures are gen- 
erated by the proliferation and eventual ordering of one di- 
mensional (ID) topological defects. A prominent example is 
provided by the stripe phase proposed to exist in the high tem- 
perature superconductors. In this case, each stripe is an anti- 
phase boundary for the antiferromagnetic order parameter. In 
general, it is difficult to prove the existence and understand the 
origin of these phases due to the complexity of the underlying 
frustrated model. 

In this paper, we will consider the Ji — J2 Heisenberg 
Hamiltonian on a square lattice with an additional term that 
makes the model quasi-exactly solvable for the fully frus- 
trated point J2/J1 = 0.5. At this point, some of the exact 
ground states are valence bond crystals with deconfined frac- 
tional S = 1/2 excitations (spinons). In addition, spin-charge 
separation occurs if the system is doped with a low concen- 
tration of holes. We also show that this particular point has 
an emergent Z2 gauge symmetry fl2il and can be associated 
with ei first order quantum phase transition between two dif- 
ferent valence bond orderings. The physical manifestation of 
the emergent gauge symmetry is a divergent susceptibility for 
the creation of ID topological defects that can be identified 
with twin-boundaries of an underlying orientational ordering. 
We will see that the common origin of these exotic behaviors 
is a dynamical decoupling of the 2D magnet into ID struc- 
tures. 

We will start by considering the following 5 = 1/2 Hamil- 
tonian on a square lattice: 

if = Ji ^ Si • Sj + J2 ^ Si • Sj 

(iJ> ((iJ» 

where and ((i, j)) denote nearest neighbors and second 
nearest neighbors respectively. The index a denotes the sites 
of the dual lattice (plaquettes) and ijkl are the four spins of 
the corresponding plaquette in cyclic order Note that the pla- 
quette interaction of H is similar to the one introduced by a 
four cyclic exchange (the only difference is in the sign of the 



LA-UR-04-44 1 9 February 2, 2008 



last term) fisll . In particular, we will consider here the fully 
frustrated point J2 = Ji/2 and K — Ji/8. For this set of 
parameters and up to an irrelevant constant the Hamiltonian 
H can be rewritten as: 

The operator projects the spin state of the plaquette a onto 
the subspace with total spin = 2. 

Exact ground states of Hp. It is clear that any state having 
at least one singlet dimer per plaquette is a ground state of Hp. 
This is because Hp is positive semi-definite and the condition 
of at least one singlet in the plaquette a implies that S*" < 1 
(there is no ~ 2 component)- The same procedure gives 
rise to the Majumdar-Ghosh [l|| and the Affleck-Kennedy- 
Lieb-Tasaki (AKLT) Q exact ground states for 5 = 1/2 
and S — \ chains respectively, and to more general ideas for 
constructing solvable 2D models fisl [T6ll . 

We found two different families of states that fulfill the 
condition of having at least one singlet dimer per plaquette. 
The first family is generated by the state which is illustrated 
in Fig^. These states are simply products of local singlet 
dimers which are represented by ellipsoids. In other words, 
the singlet dimers are completely localized and there is emer- 
gent U{1) gauge symmetry associated with this localization 
1121 . Any array of dimers along a given diagonal direction, 
(1, 1) or (1, 1), can be rotated by 7r/2 as indicated by the ar- 
rows of Fig[^. It is important to note that successive rota- 
tions have to be done along the the same diagonal direction. 
The degeneracy of this family is 2^^*+^ where Nd oc \fWs 
is the number of diagonal chains and Ng is the total number 
of sites. This degeneracy can be associated with a Z2 gauge 
emergent symmetry that changes the dimerized order param- 
eter of each individual diagonal zig-zag chain. By Z2 gauge 
symmetry we mean a local symmetry transformation that acts 
on each individual zig-zag chain mapping one of the two pos- 
sible dimerized states into the other one. In other words, this 
family of ground states is formed by configurations with par- 
allel diagonal arrays of vertical or horizontal dimers. 

There is an alternative way of viewing the local Z2 emer- 
gent symmetry. We can imagine that our system has an un- 
derlying orientational ordering given by the staggered dimer 
phase illustrated in Fig^. The two possible orientations are 
horizontal and vertical. The energy cost for creating a twin- 
boundary or inter-phase between the two different orientations 
is equal to zero. Consequently, the system has a divergent sus- 
ceptibility for the creation of ID topological defects, meaning 
that a weak coupling with another field can easily stabilize a 
particular array of twin boundaries. If the system is still in- 
variant under translations, the resulting array of topological 
defects will be periodic. The ordering of ID topological de- 
fects has been proposed in the past as a possible outcome of 
competing interactions in the high temperature superconduc- 
tors 111. 

The second family (Figs^, c and d) appears when we force 
the states to have interfaces between vertical and horizontal 
configurations along the two diagonal directions (1,1) and 
(1,1). Under this condition the state is forced to create a ds- 
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feet at the intersection between the two diagonal interfaces. In 
particular, it is possible to have an 5 = 1/2 defect or localized 
spinon as is illustrated in Fig[ni. It is easy to prove that there 
is no ground state with more than one localized spinon. Since 
the position of the defect is arbitrary, the degeneracy of each 
of these configurations is proportional to Ns. 




FIG. 1: Different ground states of Hp. The ellipsoids represent sin- 
glet dimers. 

Deconfined fractional excitations. What are the low energy 
excitations of the ground states of Fig[0' Let us first consider 
the family of solutions illustrated in Fig[3i. As we can see 
in FiglJt, if we excite one singlet dimer into a triplet state, 
the two parallel 5 = 1/2 spins can propagate along the di- 
agonals without an energy cost proportional to the separation 
between them. Consequently, the two spinon excitations are 
deconfined. Note that when the spinous propagate along one 
diagonal zig-zag chain they do not "see" the other chains be- 
cause the two possible dimerized configurations (horizontal 
and vertical) have exactly the same energy. In other words, 
the effective dimensionality of the low energy spectrum is dy- 
namically reduced from D=2 to D=l. The emergent Z2 gauge 
symmetry that we mentioned above is the mathematical man- 
ifestation of this dynamical decoupling. The Z2 dimerized 
order parameter of each diagonal zig-zag chain is decoupled 
from the corresponding order parameter of the other chains. 
The most notorious physical consequence is the emergence of 
fractional excitations which are characteristic of ID systems. 

Note that the previous analysis is only valid when the 
spinous are moving along the two diagonal directions. If, 
for instance, they propagate horizontally, they will "feel" the 
confining interaction characteristic of 2D systems because 
the horizontal chains are not dynamically decoupled. Con- 
sequently, although our system is 2D, its low energy spinon 
excitations are free to move only along ID paths. The same 
type of deconfined spinon excitations are obtained for the sec- 
ond family of ground states. 
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A similar situation occurs when one hole is introduced 
(FigEb)- The charge and the spin of the added hole get decon- 
fined because there is no energy cost for the string generated in 
between the two excitations. In other words, we can expect a 
non-Fermi liquid behavior when a magnetic system described 
by Hp is doped away from half-filling. All of these "exotic" 
behaviors are just a consequence of a dynamical decoupling 
which is signaled by an emergent Z2 gauge symmetry and 
that only occurs at the point under consideration: J2 = Ji/2 
and K ~ Ji/8. To understand what can be the physical role 
of this point we need to move away from it. 



s) Two spinous 



t>) Holon and. spinon 





HH 











































































































FIG. 2: Deconfined fractional excitations, a) Two deconfined 
spinons. b) Spin-charge separation when a low concentration of 
holes is introduced in the system. 

First order quantum phase transition. The special point de- 
scribed by Hp can be easily converted into ?l first order quan- 
tum phase transition point. Note that the two configurations of 
Fig|3]are the periodic ground states with the shortest periods 
within the set shown in FigQ Therefore, they are the leading 
candidates to remain as ground states when we depart from 
the Hp point. For instance, we can add a term to the Hamil- 
tonian which favors the staggered dimer ordering shown in 
Fig|3^ when the coupling constant g is negative or the zig- 
zag dimer configuration of Figl^J) when g is positive. There 
are different realizations of such a term and we will not focus 
in any particular one. The two different dimer phases shown 
in Fig|3]break simultaneously the translation and the rotation 
symmetry of the square lattice. The first dimer configuration 
(FigEt) is four-fold degenerate while the second one (Figl^") 
has an eight-fold degeneracy. The symmetry order parameters 
for each of these phases are: 



for the staggered dimer ordering of Fig|3t and 



(3) 



Qrj 



(4) 



where Q = (7r,7r), Q = (tt, — tt), and -q — {x,y}. Note 
that {Vl\V-l^) = (±1,0) or (0,±1) for the four equivalent 
ground states of the staggered dimer phase. For the zig-zag 
phase, the non-zero component of takes the four 

possible values 1,— i} that are necessary to identify 

the eight equivalent configurations. The remaining symmetry 
group of the zig-zag phase, Q^^-, is a subgroup of tj^*, the sym- 
metry group of the staggered dimer phase. The level crossing 
that occurs at g = between the staggered and the zig-zag 
dimer states gives rise to a first order quantum phase transi- 
tion. In other words, at this point there is a discontinuous 
change of the order parameters P*'* and I?^^ (see FiglS];)- As 
we mentioned above, this level crossing is accompanied by 
the softening of the twin-boundary defects of the staggered 
dimer phase. Hence, we can think of the zig-zag dimer phase 
as a "condensation" of these twin-boundaries. 
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FIG. 3: a) Staggered dimer phase (SD). b) Zig-zag dimer phase 
(ZD), c) The two order parameters change discontinously at 3 = 
indicating a first order quantum phase transition between the phases 
illustrated in a) and b). 

What is the general feature that gives rise to these exotic 
behaviors? To answer this question it is convenient to think of 
our system as an array of diagonal zig-zag chains. For each of 
these chains we can introduce the usual Z2 dimer order param- 
eter Vi, where I is the chain index. We can now build different 
order parameters for the 2D system by choosing different peri- 
odic configurations of Di with a well defined inter-chain wave 
vector k. For instance, in the case under consideration the 
dimer phase of Fig|3t corresponds to the uniform (fc = 0) con- 
figuration Vi ~ D, i.e., the one-dimensional Z2 order param- 
eter points in the same direction for all the different chains. In 
contrast, the zig-zag dimer ordering of Fig|3j5 corresponds to 
Vi = {—lyD, i.e., there is a staggered (k ~ tt) inter-chain 
ordering. Since the different states built in this way have dif- 
ferent wave vectors k, the quantum phase transition between 
two of them is of first order At the transition point, both 
states have the same energy meaning that there are adjacent 
chains for which the two relative orientations of Di are de- 
generate. In other words, the effective coupling between these 
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two chains has been reduced to zero and a Z2 gauge symme- 
try emerges at the transition point. In general, we can say 
that these are quantum phase transitions between two broken 
symmetry states whose order parameters have something in 
common: given a particular decomposition of our 2D system 
into ID chains described by some macroscopic variable, both 
order parameters characterize different inter-chain orderings. 

For real systems, we do not expect the dynamical decou- 
pling into ID systems to be perfect. In general, there is al- 
ways some residual interaction that makes this coupling weak 
but non-zero. A similar situation occurs in the materials that 
provide an experimental realization of one dimensional sys- 
tems. The structure of these materials contains weakly cou- 
pled chains and the coupling becomes relevant when temper- 
ature is lower than some characteristic value To . However, 
the system still behaves as an array of decoupled chains for 
T > Tq. We expect the same behavior for 2D systems that get 
dynamically deconfined into ID strucutures. 

These ideas can be extended to other lattices. For instance, 
for a honeycomb lattice we can consider a Hamiltonian that is 
the sum of operators that project the the spin state of the 
plaquette a onto the subspace with total spin 5° = 2, 3. As 
in the previous case, the system gets dynamically deconfined 
into chains that have two possible dimerized states. Conse- 
quently there is an emergent Z2 gauge symmetry that gives 
rise to the exotic behaviors which are described above. 

Conclusions. In summary, we introduced a simple exten- 
sion of the Ji — J2 Heisenberg model and we obtained differ- 
ent exact ground states for the fully frustrated point J2/J1 = 
0.5, K = Ji/8. Both the ground states and their low energy 
excitations exhibit exotic behaviors like the softening of ID 



toplogical defects and the emergence of deconfined fractional 
excitations. When some of the spins are replaced by holes, the 
phenomenon of spin-charge separation occurs in the limit of 
low concentration of holes. In addition, we showed that the 
point under consideration can be easily identified with a first 
order quantum phase transition. 

The emergence of deconfined fractional excitations was re- 
cently proposed to occur at a critical point by different authors 
S lioll . In particular, our model is an isotropic version 
of the anisotropic Confederate Flag model studied by Tsvelik 
Jioll . By analyzing the four chain model, he finds that there 
are two VBS separated by an approximately (1h-1) Lorentz 
invariant quantum critical point. However, the exponent that 
he obtains for the average dimerization is quite small, indicat- 
ing that in the limit of infinite number of chains the transition 
could become first order as in the isotropic model discussed 
in this paper Our results then show that this phenomenon 
is not restricted to quantum critical points. First order quan- 
tum phase transition points can also produce strong deviations 
from the normal properties of 2D systems. In our case, the 
common origin of these deviations is a dynamical decoupling 
of the 2D magnet into ID systems. Such a dynamical decou- 
pling was proposed for 2D strongly correlated models in the 
context of the high temperature superconductors itisll . As is 
well known, this effective reduction in the dimensionaUty of 
the system gives rise to completely different properties: the 
usual Fermi liquid is replaced by a Luttinger liquid and the 
magnetic spectrum is dominated by 5 = 1/2 spinon excita- 
tions. 
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